Some non-linear s.p.d.e.'s that are second order in time 

Robert C. Dalang^ and Carl Mueller^ 

Departement de Mathematiques 
Ecole Polytechnique Federale 
1015 Lausanne, Switzerland 
E-mail: robert.dalang@epfl.cli 

Department of Mathematics 

University of Rochester 
Rochester, NY 14627 USA 
E-mail: cmlr@troi.cc.rochester.edu 



Abstract 

We extend Walsh's theory of martingale measures in order to deal with hyperbolic 
stochastic partial differential equations that are second order in time, such as the 
wave equation and the beam equation, and driven by spatially homogeneous Gaussian 
noise. For such equations, the fundamental solution can be a distribution in the sense 
of Schwartz, which appears as an integrand in the reformulation of the s.p.d.e. as a 
stochastic integral equation. Our approach provides an alternative to the Hilbert space 
integrals of Hilbert-Schmidt operators. We give several examples, including the beam 
equation and the wave equation, with nonlinear multiplicative noise terms. 



1 Introduction 

The study of stochastic partial differential equations (s.p.d.e.'s) began in earnest following 
the papers of Pardoux |T^, ||l5l, |T^, and Krylov and Rozovskii Much of the literature 
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has been concerned with the heat equation, most often driven by space-time white noise, 
and with related parabohc equations. Such equations are first order in time, and generally 
second order in the space variables. There has been much less work on s.p.d.e.'s that are 
second order in time, such as the wave equation and related hyperbolic equations. Some 
early references are Walsh ||22|, and Carmona and Nualart |^, More recent papers are 



Mueller [Q, Dalang and Frangos [g, and Millet and Sanz-Sole [pJ]] . 

For linear equations, the noise process can be considered as a random Schwartz distribu- 
tion, and therefore the theory of deterministic p.d.e.'s can in principle be used. However, this 
yields solutions in the space of Schwartz distributions, rather than in the space of function- 
valued stochastic processes. For linear s.p.d.e.'s such as the heat and wave equation driven by 
space-time white noise, this situation is satisfactory, since, in fact, there is no function- valued 
solution when the spatial dimension is greater than 1. However, since non-linear functions 
of Schwartz distributions are difficult to define (see however Oberguggenberger and Russo 



13|), it is difficult to make sense of non-linear s.p.d.e.'s driven by space-time white noise in 
dimensions greater than 1. 

A reasonable alternative to space-time white noise is Gaussian noise with some spatial 
correlation, that remains white in time. This approach has been taken by several authors, 
and the general framework is given in Da Prato and Zabczyk However, there is again a 
difference between parabolic and hyperbolic equations: while the Green's function is smooth 
for the former, for the latter it is less and less regular as the dimension increases. For 
instance, for the wave equation, the Green's function is a bounded function in dimension 
1, is an unbounded function in dimension 2, is a measure in dimension 3, and a Schwartz 
distribution in dimensions greater than 3. 

There are at least two approaches to this issue. One is to extend the theory of stochastic 
integrals with respect to martingale measures, as developed by Walsh p2|, to a more general 
class of integrands that includes distributions. This approach was taken by Dalang [Q. In the 
case of the wave equation, this yields a solution to the non-linear equation in dimensions 1, 2 
and 3. The solution is a random field, that is, it is defined for every (t, x) G R4- x R'^. Another 
approach is to consider solutions with values in a function space, typically an L^-space: for 
each fixed t G R+, the solution is an L^-function, defined for almost all x G R'^. This 



approach has been taken by Peszat and Zabczyk in [|I^ and Peszat ||T^ . In the case of the 
non-linear wave equation, this approach yields function-valued solutions in all dimensions. 
It should be noted that the notions of random field solution and function-valued solution are 
noi equivalent: see Leveque |1I0[ . 



In this paper, we develop a general approach to non-linear s.p.d.e's, with a focus on 
equations that are second order in time, such as the wave equation and the beam equation. 
This approach goes in the direction of unifying the two described above, since we begin in 



Section |^ with an extension of Walsh's martingale measure stochastic integral |£2[, in such 
a way as to integrate processes that take values in an L^-space, with an integral that takes 
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values in the same space. This extension defines stochastic integrals of the form 



/ G{s,--y)Z{s,y)M{ds,dy), 

where G (typically a Green's function) takes values in the space of Schwartz distributions, 
Z is an adapted process with values in L'^{Rf'), and M is a Gaussian martingale measure 
with spatially homogeneous covariance. 

With this extended stochastic integral, we can study non-linear forms of a wide class of 
s.p.d.e.'s, that includes the wave and beam equations in all dimensions, namely equations 
for which the p.d.e. operator is 

where k > 1 (see Section Indeed, in Section |^ we study the corresponding non-linear 
s.p.d.e.'s. We only impose the minimal assumptions on the spatial covariance of the noise, 
that are needed even for the linear form of the s. p.d.e. to have a function-valued solution. 
The non-linear coefficients must be Lipschitz and vanish at the origin. This last property 
guarantees that with an initial condition that is in L^(R'^), the solution remains in L^(R^) 
for all time. 

In Section ^, we specialize to the wave equation in a weighted L^-space, and remove the 
condition that the non-linearity vanishes at the origin. Here, the compact support property 
of the Green's function of the wave equation is used explicitly. We note that Peszat ^\ also 



uses weighted spaces, but with a weight that decays exponentially at infinity, whereas 
here, the weight has polynomial decay. 



2 Extensions of the stochastic integral 

In this section, we define the class of Gaussian noises that drive the s.p.d.e.'s that we consider, 
and give our extension of the martingale measure stochastic integral. 

Let r'(R'^+^) be the topological vector space of functions ip G C^(R'^"'"^), the space of 
infinitely differentiable functions with compact support, with the standard notion of con- 
vergence on this space (see Adams page 19). Let F be a non- negative and non-negative 
definite (therefore symmetric) tempered measure on R'^. That is, 

T{dx){>^*^){x)>Q, for all G ©(R'^), 

where (p{x) = ip{—x), "*" denotes convolution, and there exists r > such that 



Rd 



ndx) (Y^:^ < oo. (2.1) 
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We note that if r{dx) = f{x)dx, then 



T{dx){ip*ip){x) = dx dyif{x)f{x - y)if{y) 

(this was the framework considered in Dalang 0]). Let 5(R'^) denote the Schwartz space 
of rapidly decreasing C°° test functions, and for ip G 5(R'^), let J^ip denote the Fourier 
transform of (p: 

J^(p{ri) = / exp{—iri ■ x)(p{x) dx. 
According to the Bochner-Schwartz theorem (see Schwartz ||20|, Chapter VII, Theoreme 



XVII), there is a non-negative tempered measure /i on R such that T = jF/i, that is 



T{dx)ip{x)= fi{dr])J^ip{r]), for all G 5(R''). (2.2) 

Rrf Jw 

Examples, (a) Let 60 denote the Dirac functional. Then r{dx) = 6o{x) dx satisfies the 
conditions above. 

(b) Let < a < d and set /a(x) = x G R'^. Then fa = CaTfd-a (see Stein 

Chapter V §1, Lemma 2(a)), so T{dx) = fa{x) dx also satisfies the conditions above. 

Let F = {F{ip), ip G VilV^^^)) be an L^(fi, ^, P)-valued mean zero Gaussian process 
with covariance functional 



{ip,^P)^EiF{ip)Fm= / Tidx)iip*^lj)ix). 

As in Dalang and Frangos and Dalang 0, >— >■ F{ip) extends to a worthy martingale 
measure {t, A) i-h> Mt{A) (in the sense of Walsh ||22[, pages 289-290) with covariance measure 

Q{[0,t]xAxB) = {M{A),M{B))t = t[ T{dx)[ dy Uiy) Ib{x + y) 

Jw' ^Rd 

and dominating measure K = Q, such that 

F{^)= f [ ^{t,x)M{dt,dx), for all G P(R'^+^). 

Jr+ JRd 

The underlying filtration is (J-'t = J-'t \/ A/", t > 0), where 

J^t = cr{MsiA), s<t, Ae BbiR^)), 
J\f is the cr-field generated by P-nuU sets and Bb^Rf^) denotes the bounded Borel subsets of 
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Recall that a function (s, x, u) i-^ g{s, x; u) is termed elementary if it is of the form 

g{s,x;uj) = l]a,b]{s)lA{x)X{uj), 

where 0<a<b, Ae Bbi'R'^) and X is a bounded and jF„-measurable random variable. The 
(T-field on x R'^ x Q generated by elementary functions is termed the predictable cr-field. 

Fix T > 0. Let V+ denote the set of predictable functions (s, x; u) i-^ g{s, x; u) such that 
\\g\\+ < oo, where 

ll5'll+ = ^(/ ds T{dx) dy\g{s,y)g{s,x + y)\) . 



Recall ||22| that V+ is the completion of the set of elementary functions for the norm || ■ || + . 
For g G Walsh's stochastic integral 

Mf[A)= ! ds ! g{s,x)M{ds,dx) 

Jo J A 

is well-defined and is a worthy martingale measure with covariation measure 

Qg{[0,t]xAxB)= [ ds [ T{dx) f dylA{y)lB{x + y)g{s,y)g{s,x + y) 

Jo JRd JW 

and dominating measure 

Kg{[0,t]xAxB)= [ ds [ r{dx) [ dylA{y)lB{x + y)\g{x,y)g{s,x + y)\. 

Jo Jr'' JR.'i 

For a deterministic real-valued function {s,x) ^ g{s,x) and a real-valued stochastic 
process {Z{t,x), (t,x) G R+ x R"^), consider the following hypotheses (T > is fixed). 

(Gl) For < s < T, g{s, ■) G C°°(R'^), g{s, ■) is bounded uniformly in s, and J^g{s, ■) is a 
function. 

(G2) For < s < T, Z{s,-) G C^ilV^) a.s., Z{s,-) is jF^-measurable, and in addition, 
there is a compact set K C R'^ such that supp Z{s, ■) (Z K, for < s < T. Further, 
s I— i> Z{s, ■) is mean-square continuous from [0,T] into L^(R'^), that is, for s G [0,T], 

limi?(||Z(t,.)-^(s,-)lli.(R<^)) =0. 

(G3) Ig^z < oo, where 

Ig,z=rds[ d^E{\TZ{s,-m?) I ^^{dr^)\Tg{s,■){^-r^)\\ 

Jo JRd JH'^ 
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Lemma 1 Under hypotheses (Gl), (G2) and (G3), for all x G R'^, the function defined by 
(s, y; u) I— > g{s, x — y)Z{s, y; u) belongs to V+, and so 

Vg,z{x) = / g{s,x - y)Z{s,y)M{ds,dy) 
Jo Jr^ 

is well-defined as a (Walsh-) stochastic integral. Further, a.s., x ^ Vg^z{x) belongs to //^(R*^), 
and 

|2 



EilKzWUn^) =h,z- (2.3) 



Proof. Observe that \\g{-,x — is equal to 

E[ ds T{dz) dy\g{s,x-y)Z{s,y)g{s,x-y-z)Z{s,y + z)\) . 

Because g{s, ■) is bounded uniformly in s by (Gl), this expression is bounded by a constant 
times 

e([ ds f T{dz) I dy\Z{s,y)Z{s,y + z) 

\J0 jRd jRd 

= e([ ds f V{dz){\Z{sr)\*\Z{sr)\){-z) 

\Jo JRd 

By (G2), the inner integral can be taken over K — K = {z — y : z & K, y E iC}, and the 
sup-norm of the convolution is bounded by \\Z{s, ■)lli2(j^d), so this is 



< 



e(^I^ ds\\Z{s,-)\\l,^^,^T{K-K)j= ds\\Zis,-)\\l.^^.)r{K-K)<oo, 



by ( p.l|) and the fact that s i-^ E{\\Z{s, ■)lli2(j^<i)) is continuous by (G2). Therefore, Vg^z{x) 
will be well-defined provided we show that (s, y, u) i— g{s, x — y)Z{s, y; u) is predictable, or 
equivalently, that {s,y,u!) ^— ^ Z{s,y;uj) is predictable. 
For this, set t] = jT2-" and 

2"-l 

Zr.{s,y)=Y,Z{t],x) %tj>^j(s). 

j=0 

Observe that 

J2 / ds / T{dx) {\Z{t], ■)! * \Z{t], ■)|)(-x) 



2"-l 

< r2-"^E(||z(t;^,-)|lV))r(i^-K) 

j=0 

< oo. 



Therefore, Z„ G V+, since this process, which is adapted, continuous in x and left-continuous 
in s, is clearly predictable. Further, 

e( [ ds [ Tidx) ■) - ■)!) * i\Zis, •) - Z„is, -m-x) 

< I rf.E(||z(.,-)-^„(5,-)lli.(R.))r(K-K). 

The integrand converges to and is uniformly bounded over [0, T] by (G2), so this expression 
converges to as n — » oo. Therefore, Z is predictable. 

Finally, we prove ( p.3| ). Clearly, -E'(||fg,z||^2(i^d)) is equal to 

T r \2 



E \ dx ( / g{s,x - y)Z{s,y)M{ds,dy) 



Since the covariation measure of M is Q, this equals 



e( [ dx f ds f T{dz) f dyg{s,x-y)Z{s,y)g{s,x-z-y)Z{s,z + y)] . (2.4) 

KJR'i Jo ^R'* ^R'' / 

The inner integral is equal to {g{s,x — ■)Z{s,-)) * {g{s,x — ■)Z{s, ■)){—z), and since this 
function belongs to 5(R"') by (Gl) and (G2), (O) equals 



Ei dx ds fi{dri) \T{g{s, X — ■)Z{s, ■)){ri)\" 

\J-R'i Jo JR'i 

by ( p.2| ). Because the Fourier transform takes products to convolutions, 

J'{gis,x--)Z{s,-)){v)= [ di'e^^'-^:Fg{s,-){-O^Z{s,-){^-0, 

jRd 

so, by Plancherel's theorem, 

/ dx\Hg{s,x--)Z{s,.)){^)\^= f di'\:Fg{s,-){-O^Z{s,.){r^-a?- (2-5) 

The minus can be changed to plus, and using the change of variables ^ = rj + [r] fixed), 
we find that ( |2.3| ) holds. I 
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Remark 2 An alternative expression for Ig z is 

cT 



I,,z = ^ ds f,{dr^) \\g{s, ■) * {x,{-)Z{s, ■)) 11^2(1,.)) 



where Xrii^) = e^^'^ . Indeed, notice that is equal to 

de\j'gis,-m'mxni-)zis,-)m\', 

which, by Plancherel's theorem, is equal to 

(X»7(-)^(S, ■))lli2(Rd)- 

Fix I— ^> g{s,x) such that (Gl) holds. Consider the further hypotheses: 
(G4) ds sup^J^afi{dri)\J^g{s,-){^~v)\'^ < oo, 

(G5) For < s < T, Z{s,-) G L^(R'^) a.s., Z{s,-) is ^s-measurable, and s ^— Z{s,-) is 
mean-square continuous from [0,T] into L'^(R'^). 

Fix g such that (Gl) and (G4) hold. Set 

V = {Z : (G5) holds}. 

Define a norm || ■ ||g on V by 

\\Z\\l = Ig,Z- 

We observe that by (G4) and (G5) (and Plancherel's theorem), Ig z < Ig,z < oo, where 
4z = £ dsEi\\Z{s,-)\\l.^^.^) (^snp l^^fiidr^) \J^g(s, - v)\'^ ■ 

Let 

S = {Z eV: (G2) holds}. 

By Lemma |I|, Z i— > Vg z defines an isometry from {S, \\ ■ \\g) into L'^{Q x H'^,dP x dx). 
Therefore, this isometry extends to the closure of {S, || ■ \\g) in V, which we now identify. 

Lemma 3 V is contained in the closure of {S, \\ ■ \\g). 
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Proof. Fix ip G C(f (R'^) such that ^ > 0, the support of iJj is contained in the unit ball of 
R"^ and Jj^^ ipix) dx = 1. For n > 1, set 

il>n{x) — 'nf"il>{nx). 

Then ip^ — in »5(R°') and J-'ipniO = ^i'H/^), therefore \J-'ipn{-)\ is bounded by 1. 
Fix Z eV, and show that Z belongs to the completion oi S in || ■ H^,. Set 

Zn{s, x) = a;)l[_„_„]d(x) and •) = •) * V'm- 

We first show that Zn,m e ^, that is, (G2) holds for Zn,m- Clearly, Zn,m{s, •) G C~(R'^), 
2in,m{s, ■) is JT^-mcasurable by (G5), and there is a compact set K^,^ C R'^ such that 
supp Z„^m{s, •) C iC, for < s < T. Further, 

\\Zn,m(t, ■) — Zn^rn{s, ■)llL2(R<i) < \\Zn{t, ■) — Zn{s, ■)llL2(Rd) ^ ') ~ ^{^y ■)IIl2(R<2)) 

SO s I— > Zn,m(yS, ■) is meau-square continuous by (G5). Therefore, ^„,m G ^■ 
We now show that for n fixed, — ^n,m II5 — >■ as m — >• 00. Clearly, 

W^„,.= r^^/ diE{\J^Zr,{s,-)f)\l-J'^Pmm' [ l^{dr))\J^g{s,-M-v)\'- 

Jo ^Rd 

Because |1 — ^ ^ and 

4z„ < / dsEi I di\TZn{s,-)\^][sn^ I fi{dri)\J^g{s,C - n) 

Jo VJRd / V ^ ^Rd 

= < Ig,Z < 00, 

we can apply the Dominated Convergence Theorem to see that for n fixed. 



lim ||Z„ - Zr,,m\\g = lim \ Ig,Z„-Z„,m = 0. 



m— >oo m— >oo 



Therefore, Z„ belongs to the completion of £ in || • \\g. We now show that \\Z — Zn\\g —>■ as 
n ^ 00. Clearly, 

= / ^^«^(ll(^-^n)(s, •)|li2(R<i)) (^sup^^/i(ci?7) |JP'5(s,^-n) 

Because 

11-^ ~ '^n|lL2(R,d) < (||^||L2(Rci) + ||Z„||i2(Rci))^ < 4||Z|||2(Rd), 



and Ig^z < oo, the Dominated Convergence Theorem imphes that 

hm \\Z - ZnWg = 0, 

and therefore Z belongs to the completion of S in || • \\g. Lemma ^ is proved. I 

Remark 4 Lemma ^ allows us to define the stochastic integral Vg z = g ■ provided g 
satisfies (Gl) and (G4), and Z satisfies (G5). The key property of this stochastic integral is 
that 

E [hg,z\\L^nd)) = Ig,z- 



We now proceed with a further extension of this stochastic integral, by extending the 
map g y-^ Vg z to a more general class of g. 

Fix Z E V. Given a function s ^ G{s) G 5'(R'^), consider the two properties: 

(G6) For all s > 0, J-'G{s) is a function and 

/ (is sup / fi{d7]) \TG{s,^ — 7]) \^ < oo. 
Jo 5 Jw 

(G7) For all e C^{R'^), supo<,<T G'(s) *^is bounded on R'^. 
Set 

g = {s^ G{s) : (G6) and (G7) hold}, 

and 

n = {s^ G{s) : G{s) e C^iR'^) and (Gl)and (G4) hold}. 
Clearly, HcG. For G G Q, set 

IICIU = V Ig,z- 

Notice that Ig,z < Ig,z < oo by (G5) and (G6). By Remark |, the map G i— >• vg,z is an 
isometry from {Ti, \\ ■ \\z) into L'^{Q x R*^, dP x dx). Therefore, this isometry extends to the 
closure of (H, \\ ■ \\z) in Q. 

Lemma 5 Q is contained in the closure of (7i, || ■ \\z)- 
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Proof. Fix s ^— G{s) in Q. Let ipn be as in the proof of Lemma ^. Set 

Then Gn{s, ■) G C°^(R°') by [|g, Chap.VI, Thm.ll p.l66. By (G6), J^G„(s, •) = J^G(s)-J^V^„ 
is a function, and so by (G7), (Gl) holds for Gn- Because \J-'4'n\ < 1, (G4) holds for Gn 
because it holds for G by (G6). Therefore, Gn G H. 
Observe that 

lie -G^ III = Ig-g,^,z 

= f ds f diE{\rZ{s,-m\^ f ^{dr^)\rG{s,-){^-r^)\''\l-rM^-ri)? 
Jo Jn.'i JR.'i 

The last factor is bounded by 4, has limit as n — oo, and Ig,z < oo, so the Dominated 
Convergence Theorem implies that 



lim ||G-G„|U = 0. 

This proves the lemma. 



Theorem 6 Fix Z such that (G5) holds, and s ^ G{s) such that (G6) and (Gl) hold. Then 
the stochastic integral Vg,z = G ■ is well-defined, with the isometry property 



E (ll^G,z||i2(R,d)j - Ig,z- 



It is natural to use the notation 



vg,z = I ds i G{s,- - y)Z{s,y)M{ds,dy), 
Jo Jr'^ 

and we shall do this in the sequel. 

Proof of Theorem ^. The statement is an immediate consequence of Lemma |^. 



Remark 7 Fix a deterministic function ijj G L (R ) and set 

Xt = U, f j G{s,--y)Z{s,y)M{ds,dy) 

\ Jo JR'i I L2{R'«) 

It is not difficult to check that {Xt, <t <T) is a (real-valued) martingale. 
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3 Examples 

In this section, we give a class of examples to which Theorem ^ applies. Fix an integer k > 1 
and let G be the Green's function of the p.d.e. 

^ + (-1)'=A(% = 0. (3.1) 
As in Section 3, J^G(t){^) is easily computed, and one finds 

According to [Q, Theorem 11 (see also Remark 12 in that paper), the linear s. p.d.e 

^ + (-l)'=A(% = F(t,x) (3.2) 
with vanishing initial conditions has a process solution if and only if 

fT 



f ds [ /i(rfO \^G{s)\'^ < oo, 

Jo JlR<^ 



or equivalently, 

It is therefore natural to assume this condition in order to study non-linear forms of ( p.2|) . 
In order to be able to use Theorem |^, we need the following fact. 

Lemma 8 Suppose holds. Then the Green's function G of equation l \3. 1\ ) satisfies 

conditions (G6) and (G7). 



Proof. We begin with (G7). For e C^(R'^ 

[ I sir 



SO (G7) holds. 

Turning to (G6), we first show that 



(x«G,,.,r) = (^^^^— (3.4) 
12 



where 



oo 



< Gd,k{x) = — — / e V ^p{u,x)du, 
IKk) Jo 



7(-) is Euler's Gamma function and p{u, x) is the density of a A^(0, ul)— random vector (see 
|T9| , Section 5). In particular, 



and it is shown in and that 



(G,,fc,^) = ((l + |■n-^/i), (3.5) 

and the right-hand side is finite by ( |3.3| ). However, the proofs in |]T^ and ^ use monotone 
convergence, which is not apphcable in presence of the oscillating function x^- As in [0], 
because e"*' ' has rapid decrease, 

^i/l^_.|2)fc ' = (j" (^ (i/|^_.|2)fc j ' /^|) = ■) * fee,,,), r). 

Notice that Gd,k > 0, and so 

\p{t, ■) * {XiGd,k)\ < p{t, •) * Gd,k < e^Gd,k 

by formula (5.5) in [0, so we can use monotone convergence in the first equality below and 
the Dominated Convergence Theorem in the third equality below to conclude that 

/i) = linit^o ( (T^^^pyr, /i) = \iuitio{p{t, ■) * {xiGd,k)), T) 

= (limt|o(p(t, ■) * (x^G'rf.fc)), r) = {x(Gd,k, T), 
which proves (|3^) . Because G^jt > 0, 

sup {xiGd,k, r) < (Gd^fc, r) < oo (3.6) 

by (p. 51) and p.3|). The lemma is proved. I 



4 A non-linear s.p.d.e 

Let a : R — i> R be a Lipschitz function such that a(0) = 0, so that there is a constant > 
such that for u, Ui,U2 G R, 

!«(«)!< i^lwl and \a{ui) — a{u2)\ < K\ui — U2\- (4.1) 
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Examples of such functions are a{u) = u, a{u) = sin(M), or a{u) = 1 — e 
Consider the non-hnear s.p.d.e. 



u{t,x) + {-l)^A^^^u{t,x) = a{u{t,x))F{t,x), 



(4.2) 



^(0, x) = vo{x) 



d_ 
di 



u{0, x) = vo{x) 



where Vq G L^(R ) and vq ^ H (R ), the Sobolev space of distributions such that 



I ~ II 2 tif^f 



:|^{}o(e)r<oo. 



We say that a process {u{t,-), <t <T) with values in L^(R°') is a solution of ( [4. 21) if, for 
all t > 0, a.s., 



u{t, ■) = ^G{t) * vo + Git) *vo+ f [ G{t 

Jo Jnd 



y)a{u{s,y))M{ds,dy), (4.3) 



where G is the Green's function of ( p.l| ). The third term is interpreted as the stochastic 
integral from Theorem so (u(s, ■)) must be adapted and mean-square continuous from 
[0,T] into L2(R'^). 



Theorem 9 Suppose that l \3.3i ) holds. Then equation has a unique solution {u{t, ■), < 
t <T). This solution is adapted and mean-square continuous. 

Proof. We will follow a standard Picard iteration scheme. Set 

Mir) = ^Git) * Vo + G{t) * Vo. 
dt 



Notice that vo{t, ■) G ^^(R'^). Indeed 



L2(R'«) 



sm'm')\J^vo{OM (4.4) 



< II'"o||l2(R' 



and one checks similarly that \\G{t) * VoWi^cRd-^ < ||t'||j:/-fe. Further, t Uo(t,-) from [0, T] 
into L^ilV^) is continuous. Indeed, 



lim 



^G{t)*vo - ^G{s)*vo 



0, 



L2(Rd 
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as is easily seen by proceeding as in ( |4.4| ) and using dominated convergence. Similarly, 

lim||(:;(t)*{}o-G'(s)*{)o|| =0. 

For n > 0, assume now by induction that we have defined an adapted and mean-square 
continuous process {un{s, ■), < s < T) with values in L^(R'^), and define 

Un+l(t, ■) = Uo{t, ■) + Vn+l(t, ■), (4.5) 

where ^ 

Vn+i{t,-)= / G{t- s,- -y)a{vn{s,g))M{ds,dy). (4.6) 
Jo Jh'' 

We note that (a(M„(s, ■)), < s < T) is adapted and mean-square continuous, because by 

\\a{Un{s, ■)) - a{Un{t, ■))||l2(R<*) < K\\Un{s, ■) - Un{t, ■)\\L2(Rd), 

SO the stochastic integral in (|4.6| ) is well-defined by Lemma § and Theorem 
Set 

J{s) = sup / fiid^i) |^G(s, - r])\\ (4.7) 

By (PI), (^ and (|D, 

supo<s<r is bounded by some C < oo, so by Theorem ^ and 

using (PI) , 

^(lkn+i(t,-)|li2(R.)) < 2||no(t,-)|li2(R.) + 2 / dsE(||a(n„(s,-))|li2(R.))J(t-s) 

Jo 

< 2||«o(t,-)lli2(K.) + 2KC / rfsE(||«„(s,-)||i2(R.)). (4.8) 

Jo 

Therefore, M„+i(t, ■) takes its values in L'^ilV'-). By Lemma |10| below, {un+i(t, ■), < t < T) 
is mean-square continuous and this process is adapted, so the sequence («„, n G N) is 
well-defined. By Gronwall's lemma, we have in fact 

sup snp E{\\unit, ■)\\l2,^d)) < oo. 

0<t<TneN 

We now show that the sequence ('U„(t, ■), n > 0) converges. Let 

Mn{t) = E{\\Un+l{t, ■) - Unit, ■) (R'^)) ■ 

Using the Lipschitz property of «(■), ( [4.5| ) and ([4. 61), we see that 



M„(t) <KC [ dsMn^iis). 
Jo 
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Because supQ^^^rp Mq{s) < oo, Gronwall's lemma implies that 

oo 
n=0 

In particular, {un{t, ■), n G N) converges in X R^, dP X dx), uniformly in t G [0, T], to 

a limit ^(t, ■). Because each m„ is mean-square continuous and the convergence is uniform in 
t, ■), < t < T) is also mean-square continuous, and is clearly adapted. This process 
is easily seen to satisfy ( [4 .31) , and uniqueness is checked by a standard argument. I 



The following lemma was used in the proof of Theorem 



Lemma 10 Each of the processes {un{t,-), < t < T) defined in is mean-square 

continuous. 

Proof. Fix n > 0. It was shown in the proof of Theorem ^ that t UQ{t, ■) is mean-square 
continuous, so we establish this property for t H-i> Vn+i{t, ■), defined in (^4.6|) . Observe that 
for /i > 0, 



E{\\Vn+l{t + h,-) - Vn+lit, OIliaCRd) ) < 2(/i + h) 



where 

h = 

h = 
Clearly, 

while 



E 



E 



t+h p 2 

G{t + h — s, ■ — y)a{un{s, y))M{ds, dy) 



L2(Rd) 



-'R'' 



{G(t + h - s, ■ - y) - G{t - s, ■ - y))a{un{s, y))M{dsdy) 



L2(Rd)^ 



h < 



t+h 



dsE{\\un{s, ■)lli2(R,d))J(t + h-s), 



ds / di\T{a{u^{s,-))m? / ^l{d^l) 



Rd 

The squared ratio is no greater than 

4 



Rt* 



sin(t /i - s)|^ - ri\^) - sin((t - - r]\^) 
1^ 



/ sin(/;,|,^ — r/l*^' ^ ^ 



< 



c 



2\k ■ 



It follows that 1 2 converges to as /;.—!> 0, by the dominated convergence theorem, and 
Ji converges to because the integrand is bounded. This proves that t ^— Vn+iit,-) is 
mean-square right-continuous, and left-continuity is proved in the same way. I 
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5 The wave equation in weighted L^-spaces 



In the case of the wave equation (set = 1 in ( [4.2| )), we can consider a more general class 
of non-linearities a(-) than in the previous section. This is because of the compact support 
property of the Green's function of the wave equation. 

More generally, in this section, we fix T > and consider a function s ^ G{s) G iS'(R'^) 
that satisfies (G6), (G7) and, in addition, 

(G8) There is > such that for < s < T, supp G{s) C 5(0, R). 

Fix K > d and let 9 : R'^ — > R be a smooth function for which there are constants 
< c < C such that 

c(l A \x\-^) < e{x) < G{1 A \x\-^). 
The weighted L^-space is the set of measurable / : R'^ R such that ||/||l2 < oo, where 



T 

^0 



f\x)e{x)dx. 



Let Hn = {x eYi'^ : nR< \x\ <{n+ set 

\\f\\LHH^)=(^J^ fix)dx^ ' , 
and observe that there are positive constants, which we again denote c and C, such that 

oo oo 
n=0 n=0 

For a process (Z(s, ■), < s < T), consider the following hypothesis: 

(G9) For < s < T, Z{s,-) G Li a.s., Z{s,-) is jF^-measurable, and s y-^ Z{s,-) is mean- 
square continuous from [0, T] into L^. 

Set 

£q = {Z : (G9) holds, and there is C R*^ compact such 
that for < s < T, supp Z{s, ■) C K}. 

Notice that for Z G Se, Z{s,-) G L^(Il'^) because 9{-) is bounded below on ii' by a positive 
constant, and for the same reason, s ^ Z{s, ■) is mean-square continuous from [0,T] into 
L2(R'^). Therefore, 

VG,z= [ [ G{s,--y)Z{s,y)M{ds,dy) 
Jo >/r^' 

is well-defined by Theorem |^. 
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Lemma 11 For G as above and Z G Se, vg,z ^ o.s. and 

E{\\vGAli)< r ds\\Z{s,.)\\l.J{s) 



where J{s) is defined in (|4.7|) . 

Proof. We assume for simplicity that R = 1 and K is the unit ball in R'^. Set Dq = HqUHi 
and, for n > 1, set D„ = Ui/„Uif„+i and ■) = Z(s, ■)IdS-)- By then (G8), 



oo oo 



n=0 n=0 
oo 

Therefore, by Theorem ^ 

n=0 

< C rdsf^n-'<E{\\Z{s,.)\\l,^^JJ{s) 

JO .„_n 



n=0 



Jo * 

This proves the lemma. I 
For a process {Z{s, ■)) satisfying (G9), let \\Z\\q = {Igz)^^ where 

I'g,z = £ds\\Z{s,-)\\l.As). 

Because s i-^ \\Z{s, ■)||^2 is bounded, Iq^^ < oo provided (G6) holds. Therefore, \\Z\\g defines 

a norm, and by Lemma |TT], Z ^-^ vg,z from into x R'^, c/P x 9{x)dx) is continuous. 
Therefore this map extends to the closure of £q for || ■ H^, which we now identify. 

Theorem 12 Consider a function s i— » G{s) E 5'(R'^) such that (G6), (G7) and (G8) 
hold. Let ■), < s < T) he an adapted process with values in Lg that is mean-square 
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continuous from [0, T] into Lg. Then Z is in the closure of Sg for \\ ■ \\g, and so the stochastic 
integral vc,z is well-defined, and 

EiWvczWq) < Ig,z- (5-2) 

Proof. Set •) = Z{s, „]d(-). Then (Z„) satisfies (G9) and belongs to Sg. Because 
ll-^nl-s, ■)||i2 < ■)||2,2 and z ^5 ^^e dominated convergence theorem imphes that 



hm^^oo \\Z — Zn\\g = 0, so Z is in the closure of Sg for || ■ \\g, and (|5.2|) holds by Lemma |TT 



We now use this result to obtain a solution to the following stochastic wave equation: 

Q^u{t, x) - Au{t, x) = a{u{t, x))F{t, y), (5.3) 

Ou 

u{0,x) = vo{x), —{0,x) = vo{x). 

We say that a process {u(t,-), <t <T) with values in is a solution of ( p.3| ) if {u(t, ■)) 
is adapted, t u{t, ■) is mean-square continuous from [0,T] into Lg and 

u{t,-) = ^G{t)*Vo + G{t)*do+ [ [ G{t-s,--y)a{u{s,y))M{ds,dy), (5.4) 

Jo Jr.'' 

where G is the Green's function of the wave equation. In particular, J^G{s){^) = \^\~^ sin(t|^|) 
and (G6), (G7) and (G8) hold provided ( |3.3| ) holds with k = 1. Therefore, the stochastic 
integral in ( ^.4| ) is well-defined by Theorem 0. 

Theorem 13 Suppose 



vo G L'^{IV^), vq E H ^(R'^), and a(-) is a globally Lipschitz function. Then \5.d[ ) has a 
unique solution in Lg. 

Proof. The proof follows that of Theorem so we only point out the changes relative to 
the proof of that theorem. Because a(-) is globally Lipschitz, there is K > such that for 

U, Ml, U2 G R, 

!«(«)!< 7^(1 + |m|) and \a{ui) — a{u2)\ < K\ui ~ U2\. 
Using the first of these inequalities, ( [4.8|) is replaced by 

E(||«„+i(t, -nl.) < 2\\uo{t, OIIl^cr.) + 2KG I ds{l + E{\\un{s, 







Therefore Un+iit, ■) takes its values in Lg. The remainder of the proof is unchanged, except 
that II ■ ||L2(-Rd) must be replaced by || • ||^2. This proves Theorem I 
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